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Environmental  Barrier  Coating 


A  Combined  Thermodynamic/Kinetic  Modeling  Approach  to  Predict  SiC  Recession 


due  to  Si02  Scale  Volatility  under  Combustion  Environments 


ABSTRACT 


A  computational  approach,  which  targets  on  the  prediction  of  SiC  recession  caused  by 
SiCb  scale  volatility  under  combustion  environments,  was  developed  in  this  study.  In  this 
approach,  thermodynamic  calculation  was  integrated  with  a  gaseous-diffusion  model  to 
calculate  the  fluxes  of  volatile  species,  such  as  SiO(g),  Si(OH)4(g),  SiO(OH)2(g),  and 
SiO(OH)(g),  produced  by  the  reaction  of  SiC>2  scale  with  the  combustion  air.  The  resulted 
weight  loss  of  SiC  was  then  calculated  under  a  variety  of  combustion  environments.  The 
benefit  of  using  environmental  barrier  coating  (EBC)  in  the  protection  of  SiC  from 
recession  was  demonstrated  by  the  calculation.  It  is  shown  that  the  weight  loss  of  SiC- 
based  ceramics  could  be  significantly  reduced  when  EBCs,  such  as  mullite  (AESEOb  or 
written  as  SAECE^SiCE)  or  SrAS2  (SrO-AECE^SiCE),  are  used.  The  effects  of 
combustion  conditions,  such  as  temperature  and  total  pressure,  on  the  volatility  of  SiCE 
scale  were  also  discussed. 


1.  INTRODUCTION 

Ceramic  matrix  composites  (CMCs),  which  combine  reinforcing  ceramic  phases  with  a 
ceramic  matrix,  create  materials  with  superior  properties.  Si-based  ceramics,  such  as 
silicon  carbide  (SiC)  fiber-reinforced  SiC  ceramic  matrix  composites  (SiC/SiC  CMCs), 
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are  potential  candidates  for  turbine  engine  applications  owing  to  their  high  temperature 
strength  and  durability,  as  well  as  low  density.  However,  their  usage  as  turbine  engine 
hot-section  components  is  stumbled  due  to  their  lack  of  environmental  durability  in  high 
velocity  combustion  environments.  The  silica  scale  (SiCT),  which  is  responsible  for  the 
excellent  high  temperature  oxidation  resistance  in  dry  air,  reacts  with  water  vapor  and 
forms  gaseous  silicon  oxide  and  hydroxide  species  [1]  under  combustion  environments. 
The  volatilization  of  these  species  results  in  further  oxidation  of  SiC  and  recession  of  Si- 
based  ceramics.  To  minimize  the  recession  of  Si-based  ceramics,  environmental  barrier 
coatings  (EBCs)  have  been  developed  to  prevent  silica  scale  from  reacting  with  the 
combustion  gases. 

The  first  generation  of  EBCs  is  based  on  mullite  (Ab,SbOn  or  written  as  SAbOiCSiCT) 
due  to  its  low  coefficient  of  thermal  expansion  (CTE),  excellent  chemical  compatibility 
with  Si-based  ceramics,  and  good  adherence  [2-4].  One  major  issue  with  the  mullite 
coating  is  its  relatively  high  silica  activity  (~0.5)  and  the  resulted  selective  volatilization 
of  silica  by  water  vapor.  To  overcome  the  recession  of  mullite,  water  vapor  resistant, 
such  as  yttria-stabilized  zirconia  (YSZ),  was  used  as  top  coat.  However,  one  critical 
weakness  of  YSZ  is  its  large  CTE,  twice  that  of  SiC  or  mullite.  Under  long-term 
exposure  to  thermal  cycling,  the  CTE  mismatch  causes  severe  cracking  and  delamination, 
and  leads  to  premature  EBC  failure  [1],  For  the  second  generation  of  EBCs,  the  YSZ  top 
coat  was  replaced  by  the  BSAS  [(l-^BaO-.rSrO-AECE^SiOi]  family  materials.  Current 
EBCs  usually  contain  three  layers,  the  Si-bond  coat,  the  mullite  intermediate  coat,  and  the 
BSAS  top  coat.  One  of  the  research  focuses  is  to  identify  new  top  coats  that  have  higher 
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temperature  capability  and  chemical/mechanical  compatibility  with  the  mullite 
intermediate  coat. 

To  accelerate  the  development  of  advanced  coating  materials  for  SiC-based  ceramics, 
integration  of  computational  modeling  approach  with  key  experiments  is  essential.  The 
computational  approach  is  used  to  identify  the  potential  materials,  while  experiments  are 
used  to  validate  the  prediction.  It  is  our  ultimate  goal  to  develop  a  modeling  tool  that  can 
be  used  to  predict  the  phase  stability  and  compatibility  of  EBCs  with  the  SiC-based 
substrate,  the  durability  and  high  temperature  capability  of  selected  EBCs,  the  volatility 
of  the  SiC>2  scale,  and  the  recession  of  SiC-based  ceramic  materials  under  a  variety  of 
combustion  environments.  The  focus  of  the  present  paper  is  to  calculate  the  volatility  of 
the  SiCE  scale  and  the  weight  loss  of  SiC  at  different  combustion  conditions  and 
understand  the  protective  roles  of  the  two  selected  EBCs.  A  thermodynamic  database, 
which  compiles  the  Gibbs  energies  of  the  gas  phase  and  the  condensed  phases  of  the 
silica-oxides  system,  was  used  to  calculate  the  equilibrium  between  the  EBCs  and  the 
combustion  gas.  The  gaseous-diffusion  model  proposed  by  Opila  et  al.  [5,  6]  was  used  to 
calculate  the  volatility  of  the  SiCE  scale,  and  will  be  briefly  reviewed  in  Section  2.1.  The 
significant  contribution  of  this  paper  is  that  multi-component,  multi-phase  equilibrium 
calculation  is  directly  coupled  with  the  gaseous-diffusion  model  to  calculate  the 
equilibrium  partial  pressures  and  the  fluxes  of  volatile  silicon  species  given  a  combustion 
condition.  More  importantly,  the  thermodynamic  database  of  the  silica-oxides  system 
allows  us  to  have  an  insight  view  on  how  a  selected  EBC  can  effectively  protect  SiC  from 
recession  through  the  calculation  of  the  activity  of  SiC>2,  as  will  be  demonstrated  in 
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Section  3.2.  The  purpose  of  this  work  is  not  to  validate  the  gaseous-diffusion  model 
proposed  by  Opila  et  al.  [5,  6],  while  the  effects  of  combustion  condition,  such  as 
temperature,  total  pressure,  and  the  equivalent  factor,  will  be  discussed  in  Section  3.3. 
Limitations  of  the  approach  will  be  discussed  in  Section  3.4,  and  conclusions  be 
presented  in  Section  4. 


2.  THERORY  AND  APPROACH 

2.1  Volatile  Kinetics  of  Si02 

SiC  is  thermodynamically  unstable  in  an  oxidizing  environment  and  forms  an  outer  scale 
of  Si02.  As  presented  in  the  paper  by  Opila  et  al.  [5],  in  a  combustion  environment 
containing  02,  C02,  and  H20,  SiC  can  be  oxidized  by  any  or  all  of  the  following 
reactions: 


SiC  +  ^02(g)  =  Si02+C0(g)  (1) 

SiC  +  3  C02  (g)  =  Si02  +  4  CO(g)  (2) 

SiC  +  3  H20(g)  =  Si02  +  CO(g )  +  3  H2(g)  (3) 

Due  to  the  slow  growth  rate,  Si02  forms  a  protective  layer  for  the  SiC  under  oxidizing 
environment.  However,  under  a  reducing  environment,  or  in  an  oxidizing/reducing  gas 
mixture,  such  as  H20/H2  or  C02/C0  mixtures,  the  Si02  can  be  reduced  to  form  volatile 
SiO(g)  by  the  following  reactions: 
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Si02+H2(g)  =  Si0(g)  +  H20(g ) 


(4) 


Si02  +  CO(g)  =  SiO(g)  +  C02(g )  (5) 

Similarly,  in  water-vapor  containing  environments,  SiCE  scale  may  react  with  FECEg)  to 
form  volatile  hydroxides  or  oxyhydroxides  by  one  of  the  following  reactions  [5]: 


Si02  +  H20(g )  =  SiO{OH)2{g)  (6) 

Si02  +  2  H20{g)  =  Si(OH)4(g )  (7) 

2 Si02  +  3 H20(g)  -  Si20(0H)6(g)  (8) 

Si02  +^H20(g)  -  SiO(OH)(g)  +  ^02(g)  (9) 

2Si02+3H20(g)  =  Si2(OH\(g)  +  ^02(g)  (10) 


It  is  seen  from  Equations  (4)-(10)  that  SiCE  volatility  is  attributed  primarily  to  the 
formation  of  gaseous  species,  such  as  SiO(g),  Si(OH)4(g),  Si(OH)x(g),  and  SiO(OH)x(g). 
Notice  that,  under  a  particular  combustion  environment,  only  a  few  of  the  above  chemical 
reactions  are  predominant.  The  partial  pressure  of  each  gas  species  can  be  calculated 
when  the  system  reaches  equilibrium.  Taking  the  reaction  shown  in  Eq.  (7)  as  an 
example,  we  have: 


P rmiii  ,  1  “A  Gf 

,  f  i  8)  =  exp( - ?-) 

pi  1  T?rT' 

rH20(g)aSi02  K1 


(11) 
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Where  A G°f  is  the  Gibbs  energy  change  of  the  chemical  reaction,  R  the  gas  constant,  and 


T  the  temperature  in  Kelvin;  Pshoh >4(g) ,  and  PH,o(g)  are  the  partial  pressures  of 
Si(OH)4(g)  and  H20(g),  respectively,  and  aSi0  the  activity  of  SiCE  in  the  scale.  It  should 

be  pointed  out  that,  the  equilibrium  partial  pressures  of  these  volatile  species,  such  as 
Si(OH)4(g)  in  the  above  equation,  are  directly  related  to  the  activity  of  SiCK  Apparently, 
low  aSi0n  will  suppress  the  chemical  reactions  going  to  the  right  hand  side,  reduce  the 

equilibrium  partial  pressures  of  volatile  species,  and  consequently  limit  the  volatilization 
of  SiOo.  This  point  will  be  elaborated  in  Section  3:  RESULTS  and  DISCUSSIONS. 

While  thermodynamics  can  only  provide  equilibrium  partial  pressures  of  these  gas 
species,  the  gaseous-diffusion  model  developed  by  Opila  and  colleagues  [5,  7]  is  used  to 
calculate  their  fluxes  in  this  study.  In  their  model,  the  fluxes  of  volatile  silicon  oxides, 
hydroxides,  or  oxyhydroxides  are  controlled  by  diffusion  through  a  gaseous  boundary 
layer.  In  laminar-flow  conditions,  the  following  equation  [8]  is  used  to  calculate  the 
boundary-layer-limited  flux  of  species  i: 


i  i 

=0.664^2  5^ 


Dpyt 

L 


(12) 


where  /,  is  the  mass  flux  of  volatile  species  i,  such  as  SiO(g)  and  Si(OH)4(g),  Re  the 
Reynolds  number,  Sc  the  Schmidt  number,  D  the  interdiffusion  coefficient  of  the  volatile 
species  in  the  boundary-layer  combustion  gas,  pv  the  density  of  the  volatile  species  i  at 
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the  solid/gas  interface,  and  L  the  characteristic  length  of  the  test  specimen  parallel  to  the 
direction  of  the  gas  flow.  Expanding  the  Reynolds  and  Schmidt  numbers  results  in, 


Ji  =  0.664(^)2  (-^-)3 

rj  pD  L 


(13) 


where  p  is  the  density  of  the  boundary-layer  combustion  gas,  v  the  linear  gas  velocity, 
and  ij  the  gas  viscosity.  In  this  equation,  L  and  v  are  experimental  inputs,  //  can  be 
found  in  the  literature  [9],  and  D  can  be  calculated  by  the  Chapman-Enskog  equation 
[10].  The  gas  densities,  p  and  pv  are  determined  by  the  partial  pressures  of  the  gaseous 

species.  Given  the  combustion  gas  environment,  temperature  and  pressure,  the  partial 
pressures  of  the  combustion  gas  species:  02,  H2,  H20,  CO,  C02,  N2,  as  well  as  those 
volatile  species  at  the  solid/gas  interface:  SiO(g),  Si(OH)4(g),  etc.,  can  be  calculated 
using  the  thermodynamic  database  for  the  gas  species.  The  gas  densities  p  and  pv  can 

then  be  calculated  from  the  equation  of  state  for  an  ideal  gas: 

( i  =  SiO(g),Si(OH)4,etc.)  (14) 

( j  =  02(g),H20,N2,etc .)  (15) 


Pv, 


m, 

RT 


stpjmj 
p  =  L~ 


RT 


where  /]  and  P;  represent  the  partial  pressures  of  the  volatile  silicon  species  and  the  gas 
species  in  the  combustion  environment,  respectively,  and  M ;  are  the  corresponding 


atomic  weights,  R  is  the  gas  constant  and  T  temperature.  By  using  the  partial  pressures  of 


gaseous  species  obtained  from  thermodynamic  calculation  and  Eqs.  (13)-(15),  the  fluxes 
of  the  volatile  silicon  species  can  be  calculated.  These  values  can  then  be  used  to 
calculate  the  weight  loss  of  SiC. 

2.2  Thermodynamic  Database  for  the  Gas  Species  and  the  EEC  Systems 

In  order  to  calculate  the  partial  pressure  of  each  gas  species  and  the  fluxes  of  volatile 
species,  thermodynamic  database  for  the  gas  species  must  be  developed.  The  gas  species 
considered  in  this  study  are  those  in  the  combustion  environment,  such  as  O2,  Ho,  CO, 
CO2,  H2O,  and  N2,  and  the  products  due  to  the  reaction  of  Si02  with  the  combustion 
environments,  such  as  SiO(g),  Si(OH)4(g),  SiO(OH)(g),  SiO(OH)2(g).  Gibbs  energy 
functions  for  these  species  were  developed  based  on  the  published  thermochemical  data 
in  the  literatures  [5,  11].  The  Gibbs  energy  of  the  gas  phase  is  described  by  the  following 
equation: 


G  =  2>’sG:  +  KT^y,  In  (yjsP)  =  X>’JG;  +  KTMy,r,)  +  ^lnP]  (16) 

5=1  5=1  5=1 

where  P  is  the  total  external  pressure  in  bar,  G°  is  the  Gibbs  energy  of  the  species  s  at  the 
reference  pressure  of  one  bar,  n  is  the  number  of  species  and  ys  is  the  molar  fraction  of 
the  species  5  in  the  gas  phase.  The  fugacity  coefficient,  ys ,  is  used  to  describe  the  non¬ 
ideal  gas,  which  is  unity  for  an  ideal  gas. 
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As  has  been  pointed  out  in  Section  2.1,  the  equilibrium  partial  pressures  of  volatile 
species  in  the  gas  phase  are  directly  related  to  the  activity  of  SiCE  in  the  scale.  The 
activity  of  SiCE  is  unity  when  SiC  is  directly  exposed  to  the  combustion  environment  and 
SiCE  scale  builds  up  on  the  surface.  This  is  not  true  when  an  environment  barrier  coating 
(EBC)  is  used,  in  which  case  the  activity  of  SiCE  should  be  calculated  using  the 
thermodynamic  database  of  the  EBC  system.  In  this  study,  thermodynamic  database  for 
the  SrO-A^CE-SiOo  pseudo-ternary  system  developed  by  Zhang  et  al.  [12]  was  used  to 
calculate  the  activity  of  SiC>2  for  mullite  (AESEO13  or  written  as  3AECE-2SiCE)  or  SrAS2 
(SrOAECE-2SiCE)  coatings.  Details  on  the  thermodynamic  models  used  to  describe  the 
phases  in  the  SrO-AECE-SiCE  system,  as  well  as  database  development,  will  be  published 
in  a  separate  paper  [12].  It  should  be  pointed  out  that  SrO-AECE-SiCE,  which  was 
refereed  to  as  SAS  by  Lee  et  al.  [1],  is  abbreviated  as  SrAS2-  This  is  to  distinguish  it  from 
the  other  two  ternary  compounds  in  this  pseudo-ternary  system  as  will  be  discussed  in 
Section  3.2. 

2.3  Software  Modules  for  Calculating  Complicated  Phase  Equilibria,  Gas-Solid 
Reaction  and  Volatilization  of  SiCE 

Software  package,  an  essential  part  of  this  study,  is  used  to  deal  with  complicated  phase 
equilibria,  gas-solid  reaction  and  SiCE  volatilization.  In  this  study,  all  the  functions  are 
built  upon  Pandat  [13],  a  software  package  for  multi-component  phase  diagram 
calculations.  The  unique  feature  of  Pandat  is  its  ability  to  automatically  find  the  stable 
phase  equilibrium  in  a  multi-component,  multi-phase  system.  The  gas-solid  reaction 
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module  is  used  to  find  the  equilibrium  between  the  gas  phase  and  the  EBC  system.  Figure 
1  is  a  flow  chart  showing  the  connection  between  different  modules  of  the  software  and 
the  thermodynamic  databases.  Thermodynamic  database  for  the  EBC  system  provides  the 
software  with  Gibbs  energy  functions  of  the  condensed  phases  to  find  the  phase  equilibria 
in  the  EBC  system.  The  calculated  phase  stability  and  activity  data,  as  well  as  the  Gibbs 
energies  of  gas  species,  are  inputs  for  the  gas-solid  reaction  module.  The  equilibrium 
partial  pressures  of  gas  species  calculated  from  the  gas-solid  reaction  module  are  then 
used  by  the  kinetic  module  for  calculating  the  fluxes  of  the  volatile  silicon  oxide, 
hydroxide,  and  oxyhydroxide  species. 

3.  RESULTS  AND  DISCUSSIONS 


3.1  SiC  Weight  Loss  When  Directly  Exposed  to  the  Combustion  Environments 

Robinson  and  colleagues  [6]  carried  out  experimental  study  on  the  SiC  recession  caused 
by  Si02  scale  volatility  under  combustion  conditions.  In  their  study,  SiC  coupons  were 
directly  exposed  to  the  combustion  environments.  The  Si02  was  first  formed  on  the 
sample  surface,  and  then  volatilized  due  to  its  reaction  with  combustion  air  to  form 
volatile  species,  such  as  SiO,  Si(OH)4,  SiO(OH)2,  SiO(OH).  The  process  follows 
paralinear  kinetics,  i.e.  a  simultaneous  parabolic  scale  growth,  concurrent  with  linear 
volatilization  of  the  scale.  However,  steady  state  was  quickly  established,  and  linear 
weight  loss  of  SiC  was  measured  in  their  study.  Table  I  lists  the  experimental  conditions 
of  their  experiments.  In  this  table,  T,  P  and  Vg  represent  the  temperature,  total  pressure, 
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and  gas  velocity  of  the  combustion  air.  The  equivalent  ratio,  represented  by  <j> ,  is  defined 


as  the  fuel-to-air  ratio  with  the  total  hydrocarbon  (fuel)  content  normalized  to  the  amount 
of  oxygen.  At  ^  =  1 ,  combustion  results  in  complete  consumption  of  fuel  and  oxygen.  For 
the  case  of  fuel  lean,  i.e.,  </>  <  1 ,  the  combustion  air  will  contain  N2,  H20,  C02,  and  02  as 
major  species;  while  for  the  case  of  fuel  rich,  i.e.,  <j)  >  1 ,  the  combustion  air  will  contain 
N2,  H20,  C02,  H2  and  CO. 


The  calculation  condition  used  in  this  work  is  listed  in  Table  II,  and  the  calculated  SiC 
loss  is  compared  with  the  experimental  data  of  Robinson’s  work  [6].  Figure  2  compares 
the  calculated  and  experimentally  measured  SiC  weight  loss  rate  under  fuel  lean 
condition  (refer  to  Tables  I  and  II).  It  shows  a  very  good  agreement  at  low  temperature 
range  (below  1500°C),  while  the  calculated  weight  loss  tends  to  be  higher  than  that 
determined  by  Robinson  etc. [6]  at  higher  temperature  (above  1500°C).  Figure  3  shows  a 
comparison  between  the  calculated  and  experimentally  measured  SiC  weight  loss  rate 
under  fuel  rich  condition  (refer  to  Tables  I  and  II),  and  less  satisfactory  results  are 
obtained.  The  calculated  weight  loss  rates  are  lower  than  those  determined  by 
experiments,  which  is  probably  due  to  the  formation  of  other  volatile  gas  species,  such  as 
Si(OH)6.  Due  to  the  lack  of  thermochemical  data  for  this  species,  its  Gibbs  energy  is 
unknown  and  therefore  not  considered  in  the  current  thermodynamic  database.  Figure  4, 
which  compares  the  calculated  and  experimentally  measured  SiC  weight  loss  at  two 
pressures,  indicates  that  the  model  prediction  works  well  at  different  pressures  under  fuel 
lean  condition.  In  fact,  the  calculation  for  the  higher  pressure  (lOatm)  agrees  with  the 
experimental  measurements  better,  even  for  high  temperature.  The  effects  of  temperature 
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and  pressure  on  the  weight  loss  of  SiC  are  shown  in  Figure  5  and  Figure  6.  It  is  apparent 
that  higher  temperature  and/or  higher  pressure  of  the  combustion  air  lead  to  heavier  SiC 
weight  loss. 

3.2  SiC  Weight  Loss  When  Mullite  or  SrAS2  Coating  is  used 

As  has  been  pointed  out,  the  fluxes  of  volatile  species  and  the  weight  loss  of  SiC  are 
directly  related  to  the  activity  of  S i O 2 -  It  is  seen  from  equation  (11)  that  low  aSl(X 

suppresses  the  chemical  reactions  going  to  the  right  hand  side  (equation  (7)),  and 
consequently  reduces  the  volatility  of  SiCF.  For  the  case  discussed  in  Section  3.1,  bare 
SiC  substrate  directly  exposed  to  the  combustion  environment  was  first  oxidized  to  form 
SiCF  scale,  therefore  aSia  =1  in  the  scale.  In  this  section,  we  will  demonstrate  that  SiC 

recession  rate  is  reduced  when  mullite  or  SrAS2  is  used  as  coating  due  to  the  lower 
activity  of  SiCF  ( aSia  <  1 )  in  the  EBC. 

Mullite  (3Al203-2Si02)  forms  in  the  Si02-Al203  pseudo-binary  system  as  shown  in 
Figure  7.  The  activity  of  Si02  is  calculated  along  this  pseudo-binary  at  1500K  and  plotted 
in  Figure  8.  As  can  be  seen,  aSl(X  is  1  in  the  Si02+Mullite  two-phase  field,  and  decreases 

from  1  to  0.506  in  the  single  mullite  field  and  keeps  constant  0.506  in  the  Mullite+AFCF 
two-phase  field.  Figure  9  shows  the  activity  of  SiCF  in  mullite  as  a  function  of 
temperature.  This  plot  does  not  make  sense  at  first  sight  since  the  activity  of  SiCF 
decreases  as  temperature  increases.  This  figure  needs  to  be  understood  in  combination 
with  Figure  7  and  Figure  8  together.  It  is  seen  from  Figure  7  that  mullite  has  a  certain 
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range  of  homogeneity  and  its  phase  boundary  is  tilted.  Figure  8  indicates  that  the  activity 
of  SiC>2  changes  dramatically  within  the  mullite  single  phase  region.  It  is  in  fact  not 
suitable  to  compare  the  activity  of  SiCF  at  different  temperatures  at  the  stoichiometric 
composition  of  mullite  (AlfySFOn)  since  this  composition  may  locate  in  the  SiCF+Mullite 
two-phase  field,  single  mullite  phase  field,  or  Mullitc+AFCF  two-phase  field  depending 
on  the  temperature.  The  activity  showed  in  Figure  9  is  actually  the  activity  of  SiCF  in  the 
Mullite+Al2C>3  two-phase  field,  or  the  activity  of  SiCF  following  the  phase  boundary  of 
mullite  in  equilibrium  with  AFO3,  which  is  tilted  away  from  SiCF  at  higher  temperature. 
This  is  why  the  activity  of  SiCF  decreases  as  the  temperature  increases. 

Figure  10  is  the  isothermal  section  of  the  SrCFAFCVSiCF  pseudo-ternary  at  1350°C.  The 
activity  of  SiCF  with  SrAS2  (SrO-A^CF^SiCF)  composition  is  calculated  as  a  function  of 
temperature  as  shown  in  Figure  11.  It  is  seen  from  Figure  9  and  Figure  1 1  that  the  activity 
of  SiCF  varies  from  0.53  to  0.43  in  mullite,  and  from  0.04  to  0.15  in  SrAS2  in  the 
temperature  range  of  1400K  to  1800K.  SrAS2  is  therefore  expected  to  better  protect  SiCF 
from  volatilization  and  SiC  recession  due  to  the  lower  activity  of  SiCF  in  it.  This  point  is 
clearly  shown  in  Figure  12,  in  which  the  SiC  weight  loss  rate  is  calculated  as  a  function 
of  temperature  under  the  combustion  condition  of  fuel  lean  at  6  atm.  These  calculations 
indicate  that  SiC  suffers  significant  weight  loss  when  it  is  directly  exposed  to  the 
combustion  environment,  mullite  coating  helps  to  reduce  the  SiC  weight  loss  by  about 
50%,  and  SrAS2  coating  helps  further  preventing  SiC  from  recession. 
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In  the  paper  of  Lee  et  al.  [1],  experimental  measurements  were  carried  out  to  determine 
the  weight  loss  of  SiC  when  the  BSAS  family  coatings  were  used.  The  SrAS2  (refereed  to 
as  SAS  in  Lee’s  paper)  coated  SiC  samples  were  exposed  to  50%  TLO-balance  O2 
flowing  at  4.4cm/s  at  1500°C.  The  total  pressure  of  H20  and  02  was  1  atm.  The  activity 
of  Si02  for  SrAS2  (Sr0Al203-2Si02)  is  calculated  to  be  0.145  at  1500°C,  and  the 
calculated  SiC  weight  loss  as  a  function  of  time  is  shown  in  Figure  13.  The 
experimentally  determined  weight  loss  by  Lee  et  al.  [1]  were  also  plotted  on  the  figure 
for  comparison.  Even  though  the  calculation  over  predicted  the  weight  loss  by  about 
20%,  fairly  good  agreement  is  obtained. 

Mullite  and  SrAS2  were  selected  in  the  calculation  to  demonstrate  their  protective  roles  as 
EBCs  due  to  the  fact  that  these  two  materials  have  been  studied  experimentally  for  such  a 
purpose.  In  addition  to  SrAS2,  there  are  two  other  compounds:  Sr2AS  (2Sr0-Al203-Si02) 
and  Sr6A9S2  (6Sr0-9Al203-2Si02)  in  the  SrO-ALCL-SiCL  pseudo-ternary  system  as 
shown  in  Figure  10.  To  the  best  of  our  knowledge,  their  potential  applications  as  EBCs 
for  SiC  have  not  been  studied  experimentally.  It  would  be  interesting  to  see  if  they  have 
such  a  potential  from  a  thermodynamic  point  of  view.  The  activity  of  Si02  in  both  Sr2AS 
and  SreAgSo  were  calculated  as  a  function  of  temperature  and  compared  with  that  of 
SrAS2  as  shown  in  Figure  14.  As  it  is  seen,  both  Sr2AS  and  Sr6A9S?  show  lower  activities 
of  Si02,  which  indicates  that  they  should  be  more  effective  in  protecting  Si02  from 
volatilization  and  SiC  from  recession.  This  is  especially  true  for  Sr2AS,  in  which  the 
activity  of  Si02  is  several  orders  of  magnitude  lower  than  that  of  SrAS2.  However,  it 
should  also  be  pointed  out  that  this  conclusion  is  reached  only  in  view  of  their  activities. 
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To  develop  good  coating  materials,  their  chemical  compatibility  with  Si-based  ceramics, 
coefficient  of  thermal  expansion  (CTE),  and  adherence  need  to  be  addressed  as  well. 
Nevertheless,  S^AS  and  S16A9S2  identified  by  the  thermodynamic  calculations  of  this 
study  are  worthwhile  for  future  experimental  investigation. 

3.3  Effects  of  the  Equivalent  Factor,  Temperature,  Pressure  and  Activity  of  SiC>2 
on  the  Partial  Pressures  of  Gas  Species 

In  order  to  understand  the  volatilization  of  SiOi  and  the  recession  of  SiC-based  ceramics 
under  different  combustion  environments,  it  is  essential  to  know  how  the  partial  pressure 
of  each  gas  species  changes  with  the  equivalent  ratio  temperature,  pressure  and 
activity  of  SiC>2.  In  this  section,  we  will  plot  the  calculated  partial  pressures  of  gas  species 
as  a  function  of  </>  at  different  temperature,  total  pressure  and  activity  of  SiCK  Figure  15 
is  the  one  for  the  combustion  condition  at  1200°C  with  a  total  pressure  of  6atm.  Figure  15 
(a)  shows  the  distribution  of  every  gas  species  in  the  combustion  air,  while  Figure  15  (b) 
is  for  the  volatile  silicon  oxide,  hydroxides,  oxyhydroxides  due  to  the  reaction  of  SiC>2 
with  the  combustion  air.  As  has  been  discussed  in  Section  3.1,  the  equivalent  ratio  ((j))  is 
defined  as  the  fuel-to-air  ratio  with  the  total  hydrocarbon  (fuel)  content  normalized  to  the 
amount  of  oxygen.  It  is  fuel  lean  when  the  equivalent  ratio  is  less  than  one  and  fuel  rich 
when  the  equivalent  ratio  is  greater  than  one.  As  is  seen  from  Figure  15  (a),  the 
concentrations  of  N2,  1FO,  and  CO2  remain  almost  constants  at  both  fuel  lean  and  fuel 
rich  conditions.  O2  is  one  of  the  major  species  at  fuel  lean  condition,  while  it  is  almost 
completely  consumed  at  fuel  rich  condition.  On  the  other  hand,  the  existence  of  IF  and 
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CO  is  negligible  at  fuel  lean,  but  they  become  the  major  species  at  fuel  rich  condition. 
The  partial  pressures  of  O2,  Ho,  and  CO  go  through  dramatic  change  near  the  area  where 
the  equivalent  ratio  is  one.  Although  the  concentrations  of  these  species  vary  slightly  due 
to  the  change  of  temperature  and  total  pressure,  their  distributions  as  a  function  of  (j>  keep 
the  same  trend  as  that  shown  in  Figure  15  (a).  Similar  plots  for  the  distribution  of  these 
species  at  different  temperatures  and  pressures,  therefore,  will  not  be  repeated  in  this 
paper.  In  the  following  discussions,  we  will  focus  on  the  variation  of  the  volatile  species: 
SiO,  Si(OH)4,  SiO(OH)2  and  SiO(OH),  since  they  are  the  ones  that  directly  contribute  to 
the  volatility  of  Si02. 

Figure  16  (a)  shows  the  partial  pressure  change  of  the  volatile  species  when  the  total 
pressure  jumps  from  6  atm  to  12  atm.  Figure  16  (b)  is  the  enlarged  plot  of  the  top  portion 
of  Figure  16  (a).  It  is  seen  from  these  two  figures  that  the  partial  pressure  change  of  SiO 
and  SiO(OH)  is  negligible,  while  the  partial  pressure  of  Si(OH)4  increases  by  five  folds, 
and  SiO(OH)2  by  three  folds  as  the  total  pressure  increases  from  6  atm  to  12  atm.  This  is 
consistent  with  Figure  6  which  shows  that  the  total  weight  loss  of  SiC  is  proportional  to 
the  total  pressure.  Figure  17  shows  the  partial  pressure  change  of  these  species  when  the 
temperature  changes  from  1200°C  to  1500°C.  Temperature  seems  to  have  a  big  impact  on 
the  partial  pressures  of  these  species.  The  partial  pressure  of  Si(OH)4  increases  by  about  5 
folds,  and  the  other  three  species  by  two  or  more  orders  of  magnitude  when  the 
temperature  increases  from  1200°C  to  1500°C.  In  order  to  demonstrate  the  effect  of  the 
activity  of  Si02,  partial  pressures  of  these  species  are  calculated  at  1200°C  and  6  atm,  one 
at  aSiQ  =1  (directly  exposure),  and  the  other  at  aSia  =0.05  (SrAS2  coating  at  1200°C). 
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These  two  calculations  are  compared  in  Figure  18,  which  clearly  demonstrates  that  the 
partial  pressure  of  every  volatile  species  is  reduced  by  more  than  one  order  of  magnitude 
when  SrASo  coating  is  used. 

3.4  Discussions  on  the  Limitations  of  the  Approach 

The  purpose  of  this  work  is  to  integrate  thermodynamic  calculation  with  a  gaseous- 
diffusion  model  to  predict  the  volatility  of  SiC>2  as  a  function  of  combustion  conditions 
and  environmental  barrier  coating  (EBC)  materials.  The  prediction  relies  on  the  reliability 
of  the  thermodynamic  database  for  the  gas  phase  and  the  condensed  phases  of  the  EBC 
system,  and  the  suitability  of  the  gaseous-diffusion  model.  The  accuracy  of  the  Gibbs 
energy  of  every  gas  species,  especially  those  volatile  species,  determines  the  predicted 
volatility  of  SiCF  and  the  recession  of  SiC.  As  shown  in  Figure  3,  the  calculated  weight 
loss  rates  of  SiC  are  lower  than  those  measured  under  fuel  rich  condition.  One  of  the 
possible  reasons  is  that  other  volatile  species  that  are  not  considered  in  the  database  play 
a  role.  Therefore,  a  more  comprehensive  thermodynamic  database,  which  includes  the 
Gibbs  energies  of  all  possible  gas  species,  is  needed  for  more  accurate  predictions. 
Thermodynamic  database  for  the  EBC  system,  which  is  Sr0-AE03-Si02  in  this  study,  is 
essential  since  the  activity  of  Si02  directly  determines  the  partial  pressures  of  volatile 
species  and  the  volatility  of  Si02.  This  database  needs  to  be  extended  to  a  higher  order 
system  if  the  modeling  approach  developed  in  this  study  is  to  be  used  for  other  coating 
materials.  The  gaseous-diffusion  model  used  in  this  study  was  based  on  the  assumption  of 
laminar-flow  conditions,  and  was  used  to  calculate  boundary  layer  diffusion  controlled 


18 


fluxes.  The  major  issue  with  this  model  is  that  in  some  circumstances,  such  as  very  high 
gas  velocity,  laminar  flow  is  expected  to  give  way  to  turbulent  conditions,  and  the 
equations  used  in  this  study  need  to  be  corrected.  The  limitation  of  this  model  has  been 
discussed  by  Opila  et  al.  [5]  and  will  not  be  repeated  here.  Although  with  limitations,  the 
modeling  approach  developed  with  this  study  is  a  very  useful  tool  in  understanding  the 
volatility  of  S i O 2  under  a  variety  of  combustion  conditions,  and  different  coating 
conditions. 


4.  CONCLUSION 

A  concept  of  integrating  thermodynamic  calculation  with  gaseous-diffusion  model  for  the 
prediction  of  SiC  recession  caused  by  Si02  volatility  under  combustion  conditions  is 
successfully  demonstrated  in  this  study.  In  this  approach,  thermodynamic  database  for  the 
gas  phase  and  the  condensed  phase  in  the  silica-oxide  system  is  used  to  calculate  the 
equilibrium  between  the  selected  EBC  and  the  combustion  environment.  The  calculated 
equilibrium  partial  pressures  are  then  used  by  the  gaseous-diffusion  model  to  find  the 
fluxes  of  volatile  gas  species  and  the  weight  loss  of  SiC  resulted  from  the  volatility  of 
SiCE.  The  advantage  of  this  approach  is  that  thermodynamic  calculation  is  directly 
integrated  with  SiC>2  volatile  kinetics  and  the  activity  effect  of  SiCE  for  a  selected  EBC 
can  be  instantly  taken  into  account.  As  a  result,  the  protective  efficiency  of  different 
EBCs  can  be  evaluated  under  a  variety  of  combustion  conditions,  which  provides  a 
valuable  guidance  for  the  intelligent  selection  of  EBCs  for  the  SiC-based  materials. 


19 


In  this  study,  the  volatility  of  SiC>2  and  the  recession  of  SiC  under  a  variety  of  combustion 
conditions  were  calculated  and  compared  favorably  with  the  available  experimental  data. 
It  is  demonstrated  that  the  SiC-based  ceramics  suffers  significant  weight  loss  when  it  was 
exposed  to  the  combustion  air  directly,  mullite  (3Al203-2Si02)  coating  reduces  the  SiC 
weight  loss  by  about  50%,  and  SrASi  (SrOAl203-2SiC>2)  coating  further  prevents  the  SiC 
from  recession.  Two  alloys:  Sr2AS  (2Sr0-Al203-Si02)  and  Sr6AgS2  (6SrO9Al2C>3-2Si02), 
were  identified  as  potential  EBCs  for  SiC  by  our  calculations.  The  even  lower  activity  of 
Si02  in  these  two  materials  made  them  very  promising  coating  materials  for  SiC.  The 
effects  of  the  equivalent  ratio,  temperature,  and  total  pressure  on  the  volatility  of  Si02  are 
also  discussed.  Temperature  is  found  to  have  significant  effects  on  the  equilibrium  partial 
pressures  of  the  volatile  species,  and  their  fluxes.  Higher  temperature  and/or  higher 
pressure  lead  to  heavier  weight  loss  of  SiC-based  ceramics. 
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Table  I:  Experimental  conditions  of  Robinson  et  al.  [6] 


Fuel  Condition 

T(°C) 

P  (atm) 

Vg  (m/s) 

Equivalence  Ratio 

Fuel  rich 

1225-1450 

6-6.3 

18-24 

1.69-1.97 

Fuel  lean 

1200-1450 

6-6.3 

20-23 

0.87-0.94 

Fuel  lean 

1295-1420 

10 

18-19 

0.76-0.87 

Table  II:  Calculation  conditions  used  in  this  work 


Fuel 

Condition 

T(°C) 

P  (atm) 

Vg  (m/s) 

Equivalence 

Ratio 

Constants 

Fuel  rich 

1200-1500 

6 

21 

1.8 

Fuel  lean 

1200-1500 

6 

21 

0.9 

q=5. 179x10  4  (g/cm-s) 

Fuel  lean 

1200-1500 

10 

18.5 

0.8 

D=0.287  cm2/s 

24 


Figure  1:  Flow  chart  that  shows  the  connection  of  different  modules  of  the  software  and 
the  databases 
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Figure  2:  Comparison  between  the  calculated  and  measured  SiC  weight  loss  rate  under 


fuel  lean  condition  at  6  atm  (Table  II) 
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Figure  3:  Comparison  between  the  calculated  and  measured  SiC  weight  loss  rate  under 


fuel  rich  condition  at  6  atm  (Table  II) 
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Figure  4:  Comparison  between  the  calculated  and  measured  SiC  weight  loss  rate  under 
fuel  lean  condition  at  two  different  pressures  (Table  II) 
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Figure  5:  Effect  of  temperature  on  the  weight  loss  of  SiC  at  P  =  6  atm 
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Figure  6:  Effect  of  combustion  air  pressure  on  the  weight  loss  of  SiC  at  T  =  1300°C 
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Figure  7:  Calculated  SiCF-A^Ch  phase  diagram  along  with  the  experimental  data  from 
literatures  [14-23] 
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Figure  8:  Activity  of  SiCF  along  the  SiCF-AFOs  pseudo-binary  at  T  =  1500K 
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Figure  9  :  Activity  of  SiCF  as  a  function  of  temperature  along  the  m  u  1 1  i  tc/m  u  1 1  i  tc+ A 1 2O3 
phase  boundary 
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Figure  10:  Calculated  isothermal  section  for  the  SrO-AFOs-SiCF  pseudo-ternary  system 
at  1350°C  along  with  the  experimental  data  from  Dear  et  al.  [24] 
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Figure  11:  Activity  of  SiCF  in  SrASi  (SrOAFCF^SiOi)  as  a  function  of  temperature 
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Figure  12:  Calculated  weight  loss  rate  of  SiC  under  three  conditions:  no  coating,  mullite 
coating  and  SrASi  coating 
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Figure  13:  Comparison  between  the  calculated  and  measured  SiC  weight  loss  [1]  when 
SrAS2  coating  is  used 
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Figure  14:  Comparison  of  activity  of  SiC>2  in  SrASi,  SriAS,  and  S^AgSi  as  a  function  of 
temperature 
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Figure  15:  Partial  pressure  of  gas  species  at  1200°C  and  a  total  pressure  of  6  atm 
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Figure  17:  Partial  pressure  change  of  volatile  species  with  temperature  at  P=6  atm 
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Figure  18:  Effect  of  the  activity  of  SiCF  (with  or  without  coating)  on  the  partial  pressures 
of  the  volatile  species  at  T=1200°C  and  P=6  atm 
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